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Abstract 

> , 

O . We derive the nonlinear fc^-factorization for the spectrum of jets in high-mass diffractive 

^ I deep inelastic scattering as a function of three hard scales - the virtuality of the photon Q^, 
■ the transverse momentum of the jet and the saturation scale Qa- In contrast to all other 
O . hard reactions studied so far, we encounter a clash between the two definitions of the glue in 
Q . the pomeron - from the inclusive spectrum of leading quarks and the small-/5 evolution of the 
diffractive cross section. This clash casts a further shadow on customary applications of the 
Ph| familiar coUinear factorization to a pQCD analysis of diffractive deep inelastic scattering. 
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The high-mass, or leading log-^ (LL^) evolution of Diffractive Deep Inelastic Scattering 
(DDIS) is subtle and remains under scrutiny ([11121131110011113131101111111211111, for the 

discussion in connection with the planned electron- ion collider eRHIC see fl]). When cast in 
the Good-Walker-Miettinen-Pumplin form ^3] , and upon the proper renormalization of wave 
functions, the inclusive coherent forward DDIS resembles a sum of differential cross sections 
of elastic scattering of multiparton Fock states of the photon PQ I2]- As such, DDIS is an 
inherently nonlinear functional of the color dipole scattering amplitude to be contrasted to 
the linear one in inclusive DIS [21 El- A customary treatment of DDIS as DIS off an isolated 
pomeron, viewed as a hadronic state with the well defined flux in the target, an indiscriminate 
coUinear factorization analysis of DDIS final states and conclusions on glue in the pomeron 
from such an analysis (|il7j and references therein) - all must be taken with a grain of salt, 
especially in a regime dominated by multipomeron exchanges, as e.g. nuclear targets with a 
large saturation scale. Here we present an analysis of glue in the pomeron as inferred from 
the LL-| evolution of, and jets in, high-mass DDIS which strengthens these reservations. 

In the pQCD approach to DDIS, 7*p X + p' and 7*^4 X + A', high-mass DDIS starts 
with X = qqg. Our basic points derive from the comparison of nonlinear /c^-factorization 
for jets and LL-^ evolution properties of the qqg final states at very small DDIS Bjorken 
variable f3 = Q'^/{Q'^ + M^) ^ 1, where M is the invariant mass of the diffractive state. 
As far as the LL-| evolution properties of DDIS are concerned, one may insist on the linear 
color dipole representation (0, see also the subsequent [313113) similar to that for inclusive 
DIS JW^. In 1993, such a color dipole form was derived within the double leading logarithm 
approximation (DLLA) and reinterpreted as a proof of the familiar DLLA evolution of DDIS 
[2]. Evidently, DLLA is inadequate for DDIS off nuclei (and nucleons) in the most interesting 
regime of saturation. Furthermore, a utility of the emerging operational definition of the glue 
in the pomeron, especially its /c^-factorization and unitarity properties, remain unspecified. 
Indeed, the fundamental point about inclusive DIS is that one and the same unintegrated glue 
provides the linear /c^-factorization description of both the LL^ evolution and the leading 
quark spectra ( JHl and references therein). Recently, this property has been shown to extend 
to hadron- nucleus collisions, too ^3- ^^is communication we show that for high-mass 
DDIS the same property is not borne out by a rigorous nonlinear /c^-factorization for the 
final states and the evolution of the DDIS cross section. Besides the identification of this 
clash, we report a full analytic description of the spectra of leading jets in both the photon 
and pomeron fragmentation regions of DDIS as a function of three hard scales - the virtuality 
of the photon Q^, the transverse momentum of the jet and the saturation scale Qa - which 
extends earlier limited considerations [3 13 E] ■ 

The first discussion of the qqg excitation goes back to the derivation of the color dipole 
BFKL equation |3 ED]- ^^e impact parameter space the excitation amplitude equals 
^(b, qqg) = [Sqqg{h] p, r) — Sgq(b; r)]\I'qgg(p, r), where '^qqg{p, r) is the lightcone wave function 
of the qqg state, r and p are the qq and qg dipoles, Sgqg(b;p, r) and Sgg(b;r) are the S- 
matrices of the qqg and qq Fock states at the impact parameter b (we follow the notations of 
[3 1131211 1221 )■ Inclusive forward DDIS off nucleons (N) evolves as 

= ilQWqqlQQ) + {QQgWqgg " (^qqll^d) + (1) 

t=0 

For heavy nuclei (A) the integration over the momentum transfer can be carried out explicitly 
and we consider the coherent DDIS cross section per unit area in the impact parameter space: 
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Apart from the excitation of the physical qqg state, the qqg contribution to (jT}, Q describes 
also the LL-^ evolution of the qq excitation 0. For both targets, the transverse momenta of 
jets in the qqg final state add to zero: Pg + Pg + Pg = 0. Hereafter p = Pg, q = p^, the 
rapidity gap variable xjp = x/ (3 is a, fraction of nucleon's momentum carried by the pomeron, 
z is the Feynman variable - a fraction of the photon's lightcone momentum carried by the 
quark. 

We derive nonlinear fc^-factorization directly for the observed DDIS cross section without 
the separation of the ill-defined flux of pomerons which is not borne out by the pQCD treat- 
ment of DDIS One of the building blocks is the lightcone wave function of the qq state 
of the photon and qg Fock state of the quark. For transverse photons, in terms of the QED 
splitting function Pq-yi^z) = 2Nce'jaem[z'^ + (1 ~ -2)^] one has \'^gg{z,p) — p — k)P = 

|2 



Pgy{z) \tp{e^,p) - tpie^,p - k)| , where , p) = p/(^ + P ) and 



qq\ 

zil 



and -^qg{zg,p) 
2 _ 



z)g2 + m) 

(we focus on light flavors, a simple extension to longitudinal photons and heavy flavors 
will be reported elsewhere). In high-mass DDIS gluons are soft, Zg <^ 1 
does not depend on the virtuality of the parent quark, \^qg{zg,p 
AasCp p) ^ ■'/'(/^^j P " ) where the (optional) infrared parameter /x models the 

finite propagation radius of perturbative gluons |21 EIII- The technique of Refs. [IHl ED I22j 
gives the nonlinear fc^-factorization master formulas (for nuclei we cite the leading term of 
the large- Ai'c perturbation theory, where Nc is the number of colors) 
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J d'^H,id'^K2^ih, Xjp, K;i)$(b, xip, K2)H^{ki, K2, q, p) 
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where 



(k, q, p) = Vi(/i', q){[V'i(^', P + q) - P)] + [i/^jie^ p - K + q) - V,(^', p - ^)]} 

- {q^q + '^l, (5) 



H^{ki, K2, q, p) = V'^(/^^ q)[V'j(^^ p - Ki + q) - ipjie^, P - ^i)] - {q ^ q + Ki + ^2}, 

(6) 

the free-nucleon glue /(xip,k;) is related to the qq color dipole cross section by a{x,r) = 
f d'^Kf{x, — exp(zK ■ r)], ^{h, xjp, k) = S[b; cro(xip)]5'^^^(K;) + 0(b, a;ip, k) is defined in 



3 



terms of the nuclear S-matrix for the qq dipole (THJEIIES, S[b;(To] = exp[— |cro^(b)], T(b) 
is the optical thickness of the nucleus and (Jq^xsp) = J dpKf{x-ip,n). For heavy nuclei a 
useful analytic approximation is [TH] $(b,a;ip,K) ~ (^^(b, xip)/7r((5^(b, xip) + q^)^, where 



47r" 

Nc 



as((5^)GAr(a;, (5^)T(b) is the saturation scale [TH] and Gn{x,Q^) is the inte- 



grated glue of the free nucleon. For q ^ Q\ we shall often encounter 



'^n^tu \ as(q^)G'iv(a;,q2) 
-g^(b,xip) 



2-'''--'as{Q\)GN{x,Ql] 



(7) 



The jet-integrated DDIS cross section can be cast in the color dipole form, 
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d log 4 dt 
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which defines the LL-| unintegrated glue in the pomeron via a^^xjp, P, r) = / d'^Kf^ {xjp, P, n)- 
[1 — exp(iK • r)] (a somewhat odd dimension, [cr^] = [mb ■ GeV~^] is unimportant): 

X [aK{p, p + ki) - 2K{p, p + Ki + Ka) 

+4:K{p + Ki, p + Ki + ^2) - 2i^(p + Ki, p + K2) 

-f{xjp, Kl)/(XIP, p - Ki) 2K{k2, K2 + 1^1) - K{k2, P + K2] 

-2f{xjp, p)/(Xip, Ki) \K{K2, K2 + P)+ K{k2, 1^2 " - ^(«;2, P + ^1 + ^2) 

1 Ca 
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-K 



BFKL 



fnixTP, p) + ■ y" d'^Kid^K2^^f{xjp, Kl)/(xip, K2) 



X 2/r(p + Ki, p + Ki + K2) - -f^(p + /^i, p + '^2)] 

-/(xip, Ki)/(xip, p - ki)K{k2, 1^1 + 1^2) 

-f{XjP, p)f{Xjp, Ki) K{k,2, K2 + Kl) - K{k2, P + Ki + K2' 



(9) 



where K{p, q) = li/^lfi"^, p) — i/;(yU^, q) The first form of this nonlinear fc_L-factorization result 
looks as a fusion of two pomerons, described by f^xjp, /«i,2), into the third one described by 
f^{x-sp,l3,p)- In the second form we singled out the linear BFKL evolving component with 
fD{xjp,p) = cro(a;]p)[2/(a;ip,p)-/(2)(xip,p)] which describes a'^{xjp,r) = J d'^Kfn{xJp,p)[l- 
exp{iK ■ r)] j7], where f^^\x-sp,p) = ^^;^[f ® /)(a;ip,p) and ICq = CAas/2'n:^. The inter- 
pretation of foixtp, p) as a gluon density must be taken with a big grain of salt as it is not 
positive valued, see the discussion of the (anti) shadowing properties of f^^\x,p) in Ref. [7]. 
The related LL^ result for high-mass DDIS off heavy nuclei reads (here ^^'^^ = $ (8> $) 

$^(b,Xip,/?,p) = JCoi^J d'^K,id'^K,2d'^K,3 

X 2$(b, Xjp, Kl)$(b, Xjp, K2)'^'(b, Xjp,K,3)K{p - Ki - K2,P - 1^1 - K3) 
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- 2$(b, xjp, Ki)$(b, xjp, H2)^{h, a;ip, p - k,i)K{k,3, ki + K2 + ^3) 

+ $(^)(b,XlP,p) J d'^H,id'^H,2^^^\h,Xjp,l^2)K{H,i,p + + K,2) 

)<l>(2)(b,XlP, l^2)K{p - Kl,p - K2)|. (10) 

Whether the LL^ inspired operational definitions © and (fTUjl as unintegrated glue in the 
pomeron are viable or not, can only be decided upon inspecting their utility for the description 
of diffractive final states. To this end we recall a remarkable linear fc^-factorization for the 
forward qq dijets in DIS off free nucleons 



2{27rfd(TN{r ^ qq) 
dzd'^pd'^A 



fix,A)P„iz) V(e^p)-V'(^^p-A 
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Here the distribution of the experimentally measurable dijet acoplanarity momentum, A = 
Pg + Pq, is described by precisely the same unintegrated glue f{x, A) which describes the LL^ 
evolution of DIS. In diffractive qqg final states A = — q. Following Ref. |T^ we observe that 
one can cast the DDIS cross section jHI) in the Fourier form. 



{qq\a^{xjp,P,r)\qq) = dzPq^{z) J rf^prf^q /^(a;]p, /?, q) |V'(£^ p) - '^/'(^^ P - 



:i2) 



Undoing the z, p and q integrations in ()12j] . and treating p as the quark jet momentum, i.e., 
enforcing a certain unitarity interpretation on the mathematical Fourier representation (fT^ . 
one would obtain precisely the form for the diffractive dijet spectrum. However, such 
a reinterpretation does not match the dijet cross section given by Eqs. Q, ®- The source 
of this clash is that (fT^ derives from the integral form of Q upon shifts of the integration 
variables p, q, and after such shifts their meaning of the observed jet momentum is lost. Such 
shifts signal the distortions of dipoles by multipomeron exchanges which is familiar from our 
nonlinear fc^-factorization results for the single-jet and dijet spectra PHI El 122] ^ in DDIS such 
a nonlinearity and the found mismatch between the two definitions of the glue in the pomeron 
persist already for the free-nucleon target. We conclude that the DDIS dijet spectrum, given 
by the correct unitarity cut Q, is not linear /c^-factorizable in terms of the pomeron glue Q 
operationally defined by the LL-| color dipole representation for DDIS. 

In inclusive DIS, the leading quark spectrum is linear fc^-factorizable even for nuclear 
targets: 



{27iyd(TA{r ^ qq) 
dzd^pd^h 



Pg^iz) / rf^q 0(b, X, q) V(£', P) - V'(^', P - q) • (13) 
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Repeating the considerations around (O, we would obtain for DDIS dijets the representation 
(fTSj) in terms of 0^(b, Xip, /?, q), which must be compared to (jlj) integrated over the gluon 
momentum q. Upon such a comparison, we reiterate the already made point: the spectrum 
of leading quarks in DDIS off nuclei is not linear fc^-factorizable in terms of the (nonlin- 
ear) LL-i evolution-defined glue in the nuclear pomeron, the failure of the latter to describe 
the diffractive finial states shows it is not an observable with the appropriate unitarity cut 
properties. 

Now we turn to the experimentally measurable single-jet spectra. We go directly to 
the most interesting case of nuclear targets with large saturation scale, Q\{h,x-sp) S> /x^. 



5 



and illustrate our technique with hard leading quark jets, + ^ Q\{h,x-sp), a detailed 
treatment for the free-nucleon target will be reported elsewhere. This spectrum is dominated 
by the contribution from k^, < Q\{h, Xjp), and one can expand 



•jk 



(ii(\k [Cl+ ni + H2)i[ci+ Ki + K2)k 



(q + Kl + K2j 



with Bjk = Sjk — 2pjPfc(p^ + ^. Now notice that, upon the azimuthal integrations 
(if^(p,q))= J (i^K;i(i^K;2*(b,Xip,K;i)$(b,xip,K2)i^^(Ki,K2,k,p) 



,K -q + . 



(q -K 



(14) 



(15) 



where 



- +^flog.^^ + ^ 



Q\ + cf q2 



Ql 



(16) 



vanishes at q^ ^Qa ^^d, as it was anticipated, / d'^qC'^{Q\,c[^) converges at q^ ~ Q\. 
Now we notice that to DLLA the integrand of (jl3p for inclusive DIS equals 



(j)A{h,x, q)(qi5; 



ik 



L(b,x,q) 



q^(p^ + £:^)2 -4£:^(p- q) 



(p2 + £2)2 ' '-^ (p2+e2)4 

and one can try to make a contact with DDIS in which the similar role is played by 



(17) 



+ rC^(2Q^,q^) 



2(p2 + ^2^4 



;i8) 



The two integrands, (|T7jl and (fTHjl . have different dependence on (p-q). With this reservation, 
a comparison of the two suggests 



^^{h,x^,P,q)r^-^C\2Q\,q'), 



(19) 



which has a pure higher twist behavior 0^(b, Xjp, P, q) ~ 1/q^ for q^ > Q\, to be contrasted 
with (f){h,x, q) ~ 1/q^ in inclusive DIS. 

Omitting the technicalities of the derivation for minijets, p^ < Q^, we cite our principal 
result 

<5i(b,xip) 



d da^ 



dlog^dzd^hd^p [gi(b,xip)+e2 + p2]2^ 



(20) 



which interpolates between all regimes and is one of the novelties of our paper. For hard 
jets, p^ ^ <5i(b, Xjp), the expansion (fT^ amounts to vanishing intranuclear distortions. 
Indeed, in this limit the impact parameter integration gives da^ oc / d'^hQ\(h, xjp) oc A^da^, 
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while in the opposite hmit of minijets, p^+e^ <^ <5i(b, xjp), there is an obvious strong nuclear 
suppression (cf. the discussion in pT]). 

9 \ 2 
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Ada^ 



oc 



P 



(21) 



We emphasize the importance of measuring leading quark jets at fixed as a function of 
their Feynman variable z, as that would give a handle on e"^ and the width of the plateau, 
Q\(h, Xjp) + e"^, and allow an accurate determination of the saturation scale Q\. 

Now we turn to leading jets in the pomeron fragmentation region, i.e., gluons at the 
boundary of the rapidity gap. Here one integrates over the quark momenta with the result 

J d'^p J d'^Kid'^K,2^{h,Xjp,K,i)^(h,Xjp,K2)H^{Kl,K2,q,p) 

d'^Kid'^K,2^{h, Xjp, Ki)$(b, Xjp, K2) q)V'^(/^^ q) - 2V'i(Qi, q + '^l)'0^(A^^ q) 

+ i^i{Q%q + i^i)i^iiQ\,q + i^2)]-[n{{ci + Ki + K2f)-n{{Ki-K2f)], (22) 

where we used / d'^K^{h, xjp, p— k) ^ il){Q\, p) vahd for the interesting case of heavy 

nuclei with Q\ ii^. The dependence on the virtuality of the photon Q"^ is concentrated in 



qV2.2 



n{q') = J rf^pl V(e^ p + q) - ^Pie^ p) p = 2n [-i 



2a2 



log 



- 1 



\2a\/l + a/I + — a ^ ^ 

For hard gluon jets, q^ ^ Q\, the dominant contribution would come from <^ q^, 
when one can expand fi((q + Ki + ^2)^) ~ ^(q^) + 2fi'(q^)(q ■ Ki + q ■ ^2). The first term of 
this expansion gives a rise to a contribution to (|22|) of the form 




d'^K^{h,xjp,K,) %l){ii'^,q) - %jj{Q\,q + 



+ q^ 



as{QA)GN{x,Q\) 



(24) 



Making use of Eq. (j23|) . one can readily check that the contribution from the second term, 
f2'(q^)(q- K,i 2), will have a similar large-q^ asymptotics. 

For soft gluons, q^ ^ ~ Qa^ there is an apparent small-q singularity of i/?(/i^,q). 



Upon the expansion, fi((q + Ki + ^2)^) — ^^((^1 — ^2) 



Q'{kI + Av2)(q2 + 2(q- Ki 



q ■ ^^^2) + 4k;i ■ K2), one would find that all the contributions are finite at q — > 0. For 
instance, in the DIS regime of ^ Q\, the contribution from the term oc Ki ■ k;2 will be 

proportional to | / d'^KK,j^{h,xjp, n) ['j/'j(yU^, q) — 'il^i{Q\,q + k)] | . The explicit forms of 
^7(q^) and <l>(b, x, k) allow a comprehensive study of the interplay of three scales - q^,Q\ 
and - and we obtain a simple formula, which interpolates from real photoproduction to 
DIS and from soft to hard gluons: 



d 



da^ 



9 log ^ d'^hd'^q 
, 2 



7^0 



X 



Q 



-asaemejCpNc 

asjq^ + Q\)GN{x,q^ + Q\) 
as{QA)GN{x,Q\) 



(25) 
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As far as the dependence on three scales is concerned, Eq. would hold for the free-nucleon 
target too. Notice, that compared to inclusive DIS, the scaling violation in is short of one 
power of logQ^, as it was shown already in [21 El- The gluon-jet spectrum has (i) a plateau 
in the soft region, < Q\, (ii) the 1/q^ behavior for intermediate momenta Q\ < ^ 
(for free nucleons this result is known for quite a time [HIIH]) and (iii) the 1/q^ asymptotics 
for q^ > Q^. The latter regime can be viewed as the diffraction dissociation, q ^ qg, oi the 
isolated quark from the qq Fock state of the photon, in the dipole space it corresponds to 
rare qqg configurations with ^ r^, while the bulk of DDIS comes from ^ [3 E]- 
The results for nuclear mass number dependence are instructive: in hard regions (ii) and 
(iii) the impact parameter integration gives da^ oc / (Ph Q\{h, xjp) ~ A'^^^, while in the 
plateau region da^ cx / c/^b QA^ih,xjp) ~ A^/^ for < Q\ and da^ oc / c/^b ~ A^^^ for 

^ Qa- The width of the plateau has only marginal dependence on Q^, the quantity q^da^ 
will take its maximum value at q^ ~ Q\ independent of Q^. Prior to our analytic formula 
the Qa and Q dependence of the spectrum of gluons was studied numerically [11] within 
a model which corresponds to a Gaussian parameterization for </)(b, x, k). The numerical 
trends observed in ^T] are consistent with our Eq. (j25j) . 

We come to a summary. We derived the nonlinear /c^-factorization formulas, and reported 
full analytic results, for the inclusive spectrum of leading quarks and gluons in the photon 
and pomeron fragmentation regions of high- mass DDIS, respectively. We have also derived 
the nonlinear fc^-factorization representation for the unintegrated glue in the pomeron as 
defined by the LL-^ evolution of DDIS cross section in the color dipole representation. Our 
main conclusion is that, in the general case, the LL-| evolution-defined glue does not have 
the unitarity properties appropriate for a description of the jet spectra and DDIS final states. 
This mismatch entails a breaking of the familiar coUinear factorization treatment of final 
states and is especially acute for nuclear targets with large saturation scale, but it would 
become relevant also to DDIS off free-nucleon targets at very small xip, when the saturation 
effects become substantial for DIS off free nucleons too. The issue of numerical corrections to 
the DLLA for DDIS off free nucleons needs further scrutiny in view of the improved statistics 
at HERA. 

This work has been partly supported by the DFG grant 436RUS17/138/05. 
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